Deformation properties of the BCP energy density functional 
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We explore the deformation properties of the newly postulated BCP energy density functional 
(EDF). The results obtained for three isotope chains of Mg, Dy and Ra are compared to the available 
experimental data as well as to the results of the Gogny-DlS force. Results for the fission barrier 
of 240 Pu are also discussed. 



I. INTRODUCTION 

In a recent work pj we have shown that a fully mi- 
croscopic input from nuclear and neutron Equation of 
State (EOS) calculations [2| complemented by additional 
terms accounting for finite size effects, to nuclear Den- 
sity Functional Theory (DFT) can be very successfull. 
The density functional for the ground state energy E of 
a nucleus, we considered in [l|, is of the form 

E = T + E s -°- + + Ef n f + E C (1) 

where To, E s -°-, Eq are the standard expressions for ki- 
netic energy, spin-orbit term, and Coulomb energy. The 
dependence on proton and neutron densities in E°% t is 
the one obtained from the microscopic calculation [2] 
and treated in Local Density Approximation (LDA). We 
added a Finite Range (FR) term to account for a proper 
description of the nuclear surface: 

E[J[Pn,P P ] = \Y,I fd 3 rd 3 r'p t (r)v^(r~r')pAr') 

t,f 

- 5 / d 3 r Pt (r) Pt ,(r) (2) 

t,f J 

with i=proton/neutron and jt.t' the volume integral of 
Vt,t'(r). The second term in the r.h.s. of Eq (fSJ) is in- 
troduced to preserve the nuclear matter properties of the 
microscopic calculation in the bulk. For the finite range 

form factor vt,t'(r) we made a simple Gaussian ansatz, 

_ 2 / 2 

that is Vt,t' — Vt,t'e r 1 r ° . We chose a minimum of three 
open parameters: V P:P = V n ,n = Vl, V n , P = V p , n = Vu, 
and Tq. The values of these parameters together with 
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the spin-orbit strength can be found in [.lj. They were 
adjusted to reproduce the ground-state energy of some 
spherical nuclei. For the pairing part of the interaction 
we simply take the density dependent delta force stud- 
ied in [3j for effective mass equal to the bare one. With 
this so-called Barcelona-Catania-Paris (BCP) functional 
excellent results for 161 even-even spherical nuclei with 
rms values for ground state energies and charge radii, 
comparable with the most performant functionals on the 
market 0, [H, 0, 0]> were obtained. 

In this paper we continue investigating the properties 
of the BCP functional and explore how deformation prop- 
erties of nuclei are described. We will find that the per- 
formance is again excellent, as well in comparison with 
experiment as in comparison with the results of the very 
successfull Gogny D1S force Q. We have computed the 
Potential Energy Surfaces (PES) as a function of the ax- 
ially symmetric quadrupole deformation of several iso- 
topes of the Magnesium, Dysprosium and Radium species 
with the idea in mind of covering different regions and dif- 
ferent nuclear scenarios of the Nuclide Chart. To finish 
these exploratory calculations we have considered also 
the fission barrier of the heavy nucleus 240 Pu which is 
very sensitive to quadrupole deformation properties. 



II. METHODS 

As it is customary in this kind of calculations the HFB 
equation has been recast as a minimization process of 
the energy density functional where the HFB wave func- 
tion of the Bogoliubov transformation [j| is chosen to 
minimize the energy. The variational set of HFB wave 
functions is given by means of the standard Thouless 
parametrization @. The minimization process is per- 
formed by using the gradient method as it allows an easy 
and efficient implementation of constraints by using the 
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Figure 1: The Potential Energy Surfaces (PES) of some repre- 
sentative nuclei as a function of the 02 deformation parameter 
(full line, BCP1 functional; dashed one Gogny-DlS force). In 
all the cases the zero of the curves is referred to their absolute 
minimum. 



technique of Lagrange multipliers. For the present ex- 
ploratory calculations we have restricted the calculation 
to configurations preserving axial symmetry (but reflec- 
tion symmetry breaking is allowed) in order to reduce 
the computational effort. The quasiparticle operators 
are expanded in a harmonic oscillator basis written as 
the tensor product of one oscillator in the z direction 
with oscillator frequency ui z and another in the perpen- 
dicular one and characterized by uj± . Large enough basis 
(up to 26 shells in the z direction in the case of fission) 
and optimization of the harmonic oscillator lengths have 
been used to guarantee good convergence for the values 
of all physical quantities. 

In the calculations with the Gogny force the stan- 
dard approximations have been considered, namely the 
Coulomb exchange contribution to the energy has been 
replaced by the Slater approximation and the Coulomb 
pairing field has been neglected. The two body kinetic 
energy correction has been fully taken into account. For 
details pertaining the implementation of the HFB proce- 
dure and center of mass correction with the BCP func- 
tional see Ref. [lfl 



III. RESULTS 



We have performed calculations of the potential energy 
surface (PES), as a function of the mass quadrupole mo- 
ment Q20, of several isotopes of the nuclei Magnesium, 
Dysprosium and Radium. The idea is to explore some 
representative nuclei distributed all over the Nuclide 
Chart. We have compared the results of the calculations 
with the BCP energy density functional (parametriza- 
tions BCP1 and BCP2) with those obtained by using the 
D1S Gogny force which is considered here as a bench- 
mark. For the Mg isotopes we have computed the even 
mass ones between A = 20 and A = 40 (covering the 
N=8 and N=20 neutron shell closures). For the Dyspro- 
sium isotopes we have computed the ones from A = 140 
(N=74) up to A = 170 (N=104). Finally, for the Radium 
isotopes we have considered isotopes between A = 216 
(N=128) up to A = 236 (N=148). In Fig. (TJ we have 
plotted the PES (as a function of the 02 deformation pa- 



rameter, defined as 



'-) of some representa- 



tive nuclei. As the results obtained with BCP1 and BCP2 
parametrizations are almost identical, only the first ones 
are shown. The first noticeable fact is that in all the 
nuclei studied the two curves look rather similar showing 
minima, maxima and saddle points almost at the same /?2 
values. However, the relative energy of those configura- 
tions obtained with BCP1 is different from the one with 
the Gogny force. According to the number of neutrons 
we can distinguish two different regions: the first one cor- 
responds to neutron numbers greater than the mid-shell 
value, where there are differences in the prolate side be- 
tween Gogny and BCP which extend up to rather high f3 2 
values of around 0.5. Examples of this behavior are the 
nuclei 26 Mg (N=14), 32 Mg (N=20) and 144 Dy (N=78). 
The second region corresponds to neutron numbers lower 
than the mid-shell value where the spherical configura- 
tion lies, in the Gogny case, at an energy higher than 
in the BCP results and the difference increases as N in- 
creases. This means that the deformation energy Edd, 
defined as the energy difference between the spherical 
configuration and the deformed ground state, is larger 
for the Gogny force than for the BCP functional. As 
a consequence of this behavior the excitation energy of 
the oblate minimum is higher with the Gogny interac- 
tion than with the BCP functional, reaching the differ- 
ence roughly a factor of two, as in the example shown 
of 228 Ra. Examples of this behavior shown in Fig. (U) 
are 154 Dy, 164 Dy, 228 Ra and 236 Ra. Finally, for some nu- 
clei close or at shell closure like 218 Ra shown in Fig. ([I]) 
the two PES are rather similar. These systematic differ- 
ences observed between Gogny and BCP1 are probably 
a consequence of the different surface energy coefficient 
a s (the values are a s = 17.74 MeV and 17.84 MeV for 
BCP1 and BCP2 respectively @| and a s = 18.2 MeV for 
Gogny D1S Q). Another possible explanation for the 
differences in the excitation energy of prolate and oblate 
minima could be the lower level density obtained with 
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the Gogny force ( effective mass ratio of 0.7) as compared 
with the one of the BCP functional (effective mass equal 
to the physical one) that makes shell gaps stronger (see 
below). The deformation energy E^cf also depends on 
the amount of pairing correlations (see, for example Ref. 
[lol p in such a way that the stronger the pairing correla- 
tions are the smaller the value of £def ■ It turns out that 
the particle-particle correlation energy E pp = Tt(Ak), 
which is a measure of pairing correlations, is typically 30 
to 40 % stronger in the Gogny-DlS case. However, its 
effect on E^d is not strong enough as to overcome the 
other effects mentioned above that lead to an increase 
of the deformation energy of Gogny D1S with respect to 
BCP. 

Concerning the physics, we observe in Fig. ([I]) how 
shape coexistence in 26 Mg with its oblate ground state 
appears in both the BCP and D1S calculations. We also 
observe the shoulder precursor of the deformed ground 
state after angular momentum projection in 32 Mg as well 
as the prolate ground state that develops in heavier than 
32 Mg isotopes as a consequence of neutrons populating 
the fp shell (see 38 Mg). In the case of Dy we have two ex- 
amples of shape coexistence with an oblate ground state 
( 144 Dy) and a prolate one ( 154 Dy) as well as a well de- 
formed system like 164 Dy. Shape coexistence also ap- 
pears in 142 Dy (not shown) but in this case BCP pre- 
dicts a prolate g.s. whereas Gogny predicts an oblate 
one. Finally in the lower row of Fig. ([J) we have several 
Ra isotopes ranging from the spherical 218 Ra showing an 
excited superdeformed minimum, to the reasonably well 
deformed 228 Ra and 236 Ra where the excitation energies 
of the oblate minima are quite a bit higher in the case of 
the Gogny force than with the BCP functional. 

In Fig |(2]) we have displayed the magnitude of several 
physical quantities corresponding to the ground state of 
the computed Dy isotopes, using the two BCP function- 
al and the D1S force. In the lower panel of this figure 
the two neutron separation energies are plotted. The 
typical discontinuity at the semi-magic nucleus 148 Dy 
(N=82,Z=66) is clearly observed. As usual, our predic- 
tions do not compare well with experiment around that 
point because of the impact of missing correlations in 
the binding energies of transitional (not well deformed) 
systems. In well deformed systems the agreement with 
experiment is much better for all three cases consid- 
ered here. The absolute values of the binding energies 
also agree well in all the cases and, as an example, we 
can mention that in a calculation with a basis of 15 
shells the binding energies of 160 Dy are -1305.894 MeV 
for BCP1, -1305.607 MeV for BCP2, -1304.288 MeV for 
D1S whereas the experimental value is -1309.457 MeV. 
The theoretical results do not include any kind of cor- 
relation energy beyond mean-field like the rotational en- 
ergy correction that can be estimated to be 2.18 MeV 
for BCP and 3.16 MeV for D1S. In the next panel of 
Fig ((2|) a quantity related to the mean square radius, 
namely Sr = (r 2 ) 1 / 2 - ^/ijb x 1.2 x A 1 / 3 , is plotted. 
We observe that the three theoretical predictions com- 
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Figure 2: Mean field results for the selfconsistent minimum 
of all the Dy isotopes considered. In the lower panel the 
two neutron separation energy S2N is shown against the mass 
number of the isotopes along with the experimental data. In 
the next panel the quantity Sr = (r 2 } 1/2 - ^fij^ x 1.2 x A 1/3 
is depicted. In the following panel the axis ratio r\ (see text 
for details) is represented. Finally, in the upper panel the 
deformation energy E^ c { is plotted. 
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Figure 3: Same as Fig. ([2j but for the Mg isotopes and only 
the two lower panels. 



pare rather well with experimental data and surprisingly 
the radii are closer to experiment in well deformed nuclei 
which were not considered in the original fit [l[. The dif- 
ferent result obtained in 142 Dy for BCP1 and BCP2 (and 
D1S) is due to the almost degenerate oblate and prolate 
minima in this nucleus: with BCP1 the ground state is 
prolate whereas it is oblate with BCP2 and Gogny-DlS. 
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Figure 4: Same as Fig. J3} but for the Ra isotopes. 
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In the next panel the axis ratio 77 = ((z 2 )/(ie 2 )) which 

-f- ^ +2 ) (7? > 1 for 

prolate deformation, 77 < 1 for oblate deformations and 
r/ = 1 for spherical states) is plotted. As in previous 
cases, the agreement between the three theoretical re- 
sults is very good confirming what was said in discussing 
Figure (HJ about the coincidence of the position of max- 
ima and minima of the PES. Finally in the upper panel 
the deformation energy Edef is shown. The Gogny-DlS 
values for this magnitude are systematically larger by a 
few MeV than the BCP ones as discussed previously. 

The results for the selfconsistent ground state in the 
Mg isotopes is depicted in Fig. ((3J) where we have only 
plotted the curves for SW and Sr. The two additional 
quantities shown in the case of the Dy isotopes have a 
similar behavior for the Mg nuclei and are not presented 
here. Again, we observe for the S2N a reasonably good 
agreement with experiment which is of the same qual- 
ity for the three schemes. The theoretical radii also look 
rather similar in all the theoretical approaches and com- 
pare well with the scarce experimental values. 

In Fig. (@]) we show the selfconsistent results for the 
S2N of the Ra isotopes. The results for the isotopes 
220-226R, a include octupole deformation in their ground 
state with (3$ values of the order of 0.15 for all the cases 
considered. It is not the scope of the present work to 
discuss octupole deformation in detail but it is worth 
pointing out that also this deformation multipole comes 
out quite the same independent of the force. The the- 
oretical S2N results look rather similar except for the 
nucleus 218 Ra which is close to the semi-magic 214 Ra. 
As it was discussed before, the two neutron separa- 
tion energies have a somewhat erratic behavior around 
semi-magic configurations. Most surprising is the re- 
gion around A=226-228 where the theoretical predictions 
move away from experimental values. Concerning the 
radii a good agreement between the results obtained with 
BCP, Gogny and the experiment is observed. 

In Figure © we display the potential energy surface 
(PES) as a function of the quadrupole moment corre- 
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Figure 5: Fission properties for the nucleus Pu computed 
with the two parametrizations of the BCP functional as well 
as with the Gogny-DlS force. In the lower panel the potential 
energy surfaces are depicted as a function of the quadrupole 
moment. In the upper panel the octupole and hexadecapole 
moments are presented. 



sponding to the fission process of the nucleus 240 Pu com- 
puted using the two BCP functional and the Gogny 
force. As it is customary in this kind of calculations 
fill the rotational energy correction computed with the 
Yoccoz moment of inertia has been included in the PES. 
It can be concluded from this figure (lower panel) that 
the PES obtained with the BCP1 and BCP2 function- 
al are very similar and both closely follow the shape 
of the Gogny D1S PES up to the second minimum. The 
heights of the fission barriers provided by both the BCP1 
and BCP2 functionals are lower than the one calculated 
with D1S, as it is expected from the lower surface energy 
values of BCP reported above. The smaller values of the 
fission barriers for the BCP functional results go in the 
right direction as compared to the experimental estima- 
tions [H, US] but the effect of triaxiality in the first bar- 
rier in the case of the BCP functional remains to be stud- 
ied. The impact in the spontaneous fission half life <sf of 
240 Pu is also uncertain as this quantity not only depends 
on the topology of the PES but also on the collective 
inertia that has not been considered here. Preliminary 
calculations not including triaxiallity (therefore produc- 
ing too high values for the half life) yield isF = 1.2 x 10 28 s 
for BCP1, 1.1 x 10 27 s for BCP2 and 1.5 X 10 26 s for D1S. 
Therefore it seems that the larger values of the collec- 
tive inertia obtained for the BCP functionals (consistent 
with their lower pairing correlations) somehow counter- 
act their lower fission barriers providing a longer half life. 
Although the shapes involved in fission are mainly char- 
acterized by its quadrupole moment, higher multipole de- 
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Figure 6: Spherical single particle energies for protons (left 
panel) and neutrons (right panel) for the nucleus 160 Dy. The 
results for the BCP1 functional and the Gogny force (GOG) 
are given. To help identify the levels their labels alternate 
their position with respect to the line assigned to them 



formations can be important for describing some fine de- 
tails of the nuclear dynamics. In the top panel of Figure 
([5]) the expectation values of the octupole and hexade- 
capole moments are displayed as a function of Q20 and 
striking similarity between the three results is observed. 

The similarity between the results obtained with the 
BCP functionals and the Gogny force calls for a compari- 
son of the underlying single particle structure that, as it is 
well known, is essential in determining the response of the 
system to deformation. The most thorough comparison 
would involve the analysis of the single particle energy 
plots as a function of quadrupole deformation (Nilsson 
diagrams) but this is a demanding task that is deferred 
to a longer publication. We just want to mention here 
that the most relevant feature of the Nilsson diagrams, 
namely the shell gaps determining the location of min- 
ima are rather similar for both the BCP functionals and 
the Gogny force. This is so in spite of the different ef- 
fective mass ratios m*/m in the two cases, one for the 
BCP functionals and 0.7 for Gogny-DlS, that imply a 



higher level density in the former case. To illustrate the 
different level densities we have plotted in Figure iJH) for 
the nucleus 160 Dy the single particle levels at sphericity 
(Q20 = 0) for BCP1 and Gogny-DlS both for protons 
and neutrons. The higher level density of BCP's single 
particle levels is clearly observed. Results for BCP2 are 
not included, since they are very similar to the ones of 
BCP1. 

IV. CONCLUSIONS 

We have analyzed by means of a few selected examples 
the performance of the new BCP functional in what con- 
cerns quadrupole deformation. The tests performed are 
quite demanding as they include fission barriers for Pu 
or quadrupole properties over a wide range of isotopes in 
various regions of the Nuclide Chart. We have compared 
the results of our mean field calculations with experi- 
mental data whenever possible and found good agree- 
ment comparable to the one obtained for spherical magic 
or semimagic nuclei. Other quantities not directly re- 
lated to experiment like the topology of the potential 
energy surfaces have been compared to the results of a 
well performant force and with long tradition in the field, 
namely the Gogny-DlS force and the results are compa- 
rable and of the same quality. Minor differences can be 
attributed to a slightly lower surface tension for BCP 
compared to D1S and to the different pairing interac- 
tions used. From these exploratory calculations we can 
conclude that the BCP functional with the BCP1 and 
BCP2 parametrization can be used with confidence in 
the study of nuclear properties related to deformation. 
Obviously a more thorough study of the BCP functional 
with respect to deformation has to be carried out as, e.g 
the study of octupole deformation, triaxiallity, collective 
inertia, including the moment of inertia, etc. Work in 
this direction is underway. 
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